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Abstract. The singularities of the electromagnetic field are derived to include all the 
point-like multipoles representing an electric charge and current distribution. Partial 
Q^l results obtained in a previous paper [1] are completed to represent accurately all the 

' terms included in these singularities. 

^ ! 1. Introduction 

In the cases of electrostatic and magnetostatic fields of point-like dipoles, one has 
p I the well-known procedure of introducing Dirac (5-function terms for obtaining correct 

expressions of the electric and magnetic fields defined on the entire space. The 
corresponding field expressions take the following form [2]: 

SEq 4:7160 f 3^0 

where u = r /r, and 

8: B„(r) = ^(r) + ^ ^'^ ^ gg2 ^ ^(.) + (B).,,, (2) 

In these equations, by (. . .)r^o we understand an expression in which the derivatives 
are calculated supposing r 7^ 0, representing some well-known expressions of the fields. 
The expressions from equations (II]) and ([2]) are introduced in Ref. [2] as conditions of 
compatibility with the average value of the electric or magnetic field over a spherical 
domain containing all the charges or currents inside. Another procedure for introducing 
equations ([T]) and is based on an extension of the derivative didj/ (1/r) to the entire 
space [3]: 

9.94- + (3) 

A more pedagogical and suitable approach for understanding the origin of the difference 
between the electric and magnetic cases is done in Ref. [1]. Refs. [5] and [6] contain 
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generalizations of the equations ([T]) and ([2]) to the dynamic case for oscillating electric 
and magnetic dipoles. 

The objective of the present paper is to establish the singularities of the 
electromagnetic field associated to a system of electric charges and currents assimilated 
with a point-like multipolar system. These singularities are established for an arbitrary 
multipolar order in the static case. 

In section 2, the procedure of separating the 5- form singularities of the multiple 
partial derivatives of is described. In sections 3 and 4, we separate the 6- 
form singularities of the fields corresponding to the point-like equivalent multipole 
distributions, based on the multipole expansions of the fields associated to an electric 
charge or current distribution confined in a domain V. 

We point out that the formalism presented in this paper has as mathematical basis 
the properties of the irreducible tensorial representations of the proper rotations group 
[7]. For pedagogical and larger accessibility reasons, we give an explicit calculation 
based on the properties of the tensor contractions such that the procedure has a simple 
algebraic character. A counterpart of the procedure used in the present paper could be 
represented by the technique of the spherical function expansions and the cited issues 
can be a basis for such an approach. However, the use of the spherical coordinates can 
lead to the omission of the 5-form singularity contributions, as argued in the conclusions. 



2. Some delta-function identities 

The treatment of some delta-function identities in Ref. [3] can be easily generalized to 
obtain the necessary identities for higher order derivatives of 1/r. In the points different 
from the origin O, the function 1/r is a solution of the Laplace equation: 

A- = 0, r ^ . (4) 
r 

It is well-known how one can extend A(l/r) as a distribution to the entire space, as 
solution of the Poisson equation: 

A^ = -47r5(r). (5) 

Equation ([3]) represents the next step of generalization and can be continued for arbitrary 
multiple partial derivatives of 1/r i.e. for di-^ . . . (9j„ (1/r) for arbitrary n. For r 7^ 0, this 
multiple partial derivative is given by the formula 

In this equation, the coefficient C is fully symmetric with respect to ii . . .in and can be 
written as 

[-] 

= E (-l)'(2r^ -2k- 1)!!5|,,, . . . S,,,^,,,,u,,^, . . . . (7) 

k 

In the last equation, is the integer part of (3 and by Ai^i^ i^^ we understand the sum 
over all the permutations of the symbols iq giving distinct terms. From equation it 
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is obvious that we have to consider the derivatives of 1/r extended over the entire space 
as distributions (generahzed functions) which contain singular distributions (5-functions 
for example) as separate terms. 

Let be a function F{r) and suppose the existence of the integral of the product 
F{r) (f){r), with (j){r) an arbitrary smooth function (a test function from the domain of 
the distributions), on the spherical region V^, with arbitrary radius R, delimited by the 
spherical surface S/j with the center in O: 



d^x F{r) (j){r) = lim 



d?xF{r)(p{r) 



(8) 



This integral can be expressed excluding from the domain a spherical domain of 
radius e centered in O. Writing this last limit of integrals, we can interpret the function 
F{r) as a distribution defined by 



{F{r)\.^,, 4>{r) 



lim 



d^x (F(r)),^o<^(r) 



This distribution can be extended such that its support includes the point O. A new 
term 6[e — r) F{r) can be naturally introduced by the identity 

F,{r) = e{e - r) F{r, t) + e{r - e)F(r), 

associated with the extension of the integral to the entire domain P/j: 

d^xFir 



(F(r),0(r)) 



lim 



/ d^xF{r)(j){r) 



(r) 



(9) 



Moreover, we suppose the existence of the integral ([8]) for the partial derivatives of 
F. Let us consider the partial derivative diF{r, t) and the problem of extending this 
function as a distribution. The definition ([9]) becomes: 



(F(r,t), 0(r)) 



lim 

e-)>0 



+ 



<f dSuiF{r,t)(f){r)- [ d^xF{r,t)d, 



Vr\D, 



d^xdiF{r,t)4>{r) 



(10) 



where is the sphere of radius e centered in O and the Gauss theorem was employed. 
Let us apply this definition to the derivative didj{l/r) and, for simplifying the notation, 
let 



Di^,„i^{r,t) = di, 
such that we can write 



r 



lim 

£^0 



dS Ui dj - I 



/ d^x^,-^^(t){r) + I d^xidid^- 
Jve r Jt>r\v^ r 



(11) 



Since the last integral on the domain \ represents the distributions associated 
with the F- expressions for r ^ and, in the case of the electromagnetic field, they will 
be the well-known expressions of the multipole expansions, in the following, we consider 
only that part of (-Djj) containing singular distributions with point-like support i.e.. 



:i2) 
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actually, the difference 

((A,)(o) , 0) = (A„ 0) - lim / d'x {d,d,-) 0(r) 

= lim (f dS Ui dj- (j)(r) — [ d^x dj- di(f)(r] 

By -D(ro) we denote a distribution having as support the point given by the vector Vq. 
The surface integral, 

lim / dS Ui di- Sir) = — hm / dS —ViVjSir) , 

after inserting the Taylor series of the function 0(r) and since on the sphere r = e, 
becomes [3], 

lim I dS V, d^- 0(r) = - lim / d^{v) v^Vj [0(0) + ev^, {dk(p)o + . . .] (13) 
Let us introduce the angular average: 

Particularly, we have the well-known formula [8]: 

z/, ) = / 1 °' n = 2k + l, ^^^^ 

\ (n+l)!! ^{«i*2 • • • '5i„_ii„}5 n = 2k, k = 0,l,... 

Excepting the term containing 0(0), all the terms in equation ( !T4|) are proportional to 
positive powers of e and, consequently, vanish with e — 0, such that 

lim j( dSz/,9,i0(r) = -47r 0(0) = -y5.,0(O) = -^6,, (5(r), 0(r)) . 

Considering the second integral in the right-hand side of equation (ITT]) , we can write 

lim^ d\ 9,^9,0(r) = -limj^^dr J dn{i^) J^,^[(9,0)o + rz/fc(a,a,0)o + ...]= , 

such that finally 

5A^)^^^ = -y5,,5(r-) , (16) 

i.e. the delta-singularity from equation (3) of Ref. [3]. 

Let us consider the distribution Dij^. Considering only the part having O as 
support, 

((A,fc)(o), 0) =lim[/ dSiy,d,dk-(p{r)- [ d^'x (8,8^-) d,<f){r)\ . (17) 
The surface integral becomes: 

lim (f dS Ui 8j8k- 0(r) = lim / dS i^i-^r(3z/,z/fc - Sjk)4>{r) 
and, introducing the Taylor series for 0(r), 

lini / dS" Ui 8j8k- 4>{r) = in lim / -{SuiUjUk - UiSjk) [0(0) + eui {8i(j))Q + •••]) , (18) 
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and 

lim / dS Ui djdk- 4>{f) 



4:71 {SuiUjUkiyi - UiUiSjk) {dicf))^ = in (^^S^ijSki} - ^Su6jk^ {dicj))^ . (19) 



Concerning the integral on from equation ( 1171) . we have to observe that, beginning 
from this derivative order, there is a non-zero contribution for e — )■ [3]. Indeed, 
introducing equation (fT6l) in equation (|T71) and noticing that the term (9j9fc(/(r)/r))j.^o 
gives a null contribution to the limit for e — )■ 0, we can write 

-lim / d^xdjdk- di(j){r) = ^ / d'^x 6jk6{r)di(f){r) = ^ 6jk{di(f))o . 

Finally, equations ( IT71) . (ITSl) and ( IT9l) give 

{Dijk)(o) = — -S{ijdk}S{r) , (20) 

i.e. the delta-singularity from equation (4) of Ref. [3]. 

Obviously, this procedure becomes very complicated for higher order derivatives. 
Fortunately, for the electromagnetic field, some invariance properties allow a 
considerable simplification of such calculations. 

3. Singularities of the electrostatic field 

Let us consider the multipole expansions of the electrostatic field. Given an electric 
charge distribution with support included in the domain V, the scalar potential is 
expressed in the exterior of a sphere containing this domain by the following multipolar 
series: 

1 f— 11" P- • 1 (—1)"- 

Hr) = E M^^n . . . 5..^ = E M^V-II^— . (21) 

In this expansion, the coordinate system origin O is supposed in T) and P'-"^ is the n-th 
order electric multipolar moment defined by the Cartesian components in the general 
dynamic case: 

Pi,...^At) = l^d'xx,,...x,^pir,t): pW(t) = |d\rXr,t). (22) 

V 

In equation (|2TI) we employed the following notation for tensorial contractions: 

{^ji - j„_mji---j,„-Bii - im 1 n > m 
An-uBn-u ,n = m . (23) 

Ajl---jnBji---j„ii---im-n ■, U < Vfl 

For the multipole expansion of the electric field E{r) = —V $(?"), we can write 

1 f — n"-! pin) 1 f — iyra-l p. . 

Eir) = ^ E ^^V"+^||^— = -^e, E ^^9.9n • • - S.^^, (24) 
4vr£o^i n! r AttEq n\ r 

where, for simplicity, the electric charged system is considered neutral (Q = 0). We have 
to search the singularities of E given by equations (12^ . It appears that cumbersome 
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calculations are involved for higher n if we apply the formulae for higher order derivatives 
of 1/r as in the previous section. However, we can employ an invariance property 
of the electrostatic field to the substitutions of all moments P^"\ for all n, by their 
corresponding symmetric and trace-free STF projections T'^") [9l[T0l[Tl]. Retaining the 
notation p for the first order moment, this invariance stands for the invariance of the 
multipole expansion of the electrostatic field to the following substitutions: 

p, P(2), P(3), ... ^ p(2),p(3)... - 

(25) 

These STF tensors can be expressed by the following formula: 

" (2n - 1)!! Jv ^ ' *i ^ > 

which, actually, differ from the projections by numerical factors. Here, 
B^^^ . = B- B- B- B- 

Some care is necessary when one considers this invariance property when we have to 
establish the delta-type singularities of the electromagnetic field, since this property is 
true only for r 7^ 0. Indeed, as we firstly see in the electrostatic case, the substitutions 
p(") _i. -p(") in equation ( 12^ give additional terms containing A(l/r) and their 
derivatives, which, in the case r 7^ 0, can be eliminated. But, extending the expressions 
associated to the multipole expansions to the entire space, including the origin O, these 
additional terms give notable contributions. This process of searching the delta-form 
singularities of the electromagnetic field is lost in the previous paper [1] and here, we 
try to correct it. 

Let us consider the delta- singularity corresponding to the electric dipolar field: 

< = 3^e.9.9,f = -^p.W. (27) 



a result obtained by directly applying equation (1161) . 

For the 4-polar term from E, we firstly consider the expansion expressed by 
the primitive moments p("): 

E^'\r) = --^eAdjdk^ . (28) 
StiEq r 

With the help of equation (!20|) . we obtain: 

-Sjoj = Y^Si PjkS{ijdk}S{r) = j^f^i Pjk {5ijdk5{r) + 5ikdj5{r) + 5,jkdi5{r)) 

^ -ei {2P,jB,5{r) + P,jdi5{r)) = -^e^ {p^J^j5{r) + \Pj^B,5{r)] . (29) 



10£o ^ - ' ■ ' V 5£o V ' ^ ' 2 

Note that P'-^^ is symmetric. This result suggests the possibility to facilitate the 
calculation of the delta-type singularities of the electromagnetic field if instead of the 
"primitive" tensors P*^"^ we can employ the STF projections 7"^"). For the present 
case, it will be simpler to calculate the contractions of the last type of tensors with the 
coefficients C^"'* from equation (Ej). 
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For n = 2, we search the trace-free part of the symmetric tensor P(^) as 

Pi, = Vij + 6ij A . (30) 
The parameter A is determined such that Vu = 0. One obtains 

A=^Pn. (31) 



Let us introduce equation (1301) in equation ( 12811 : 

E^^\r) = —SiVikdidA- - -^e^^^ ( A- 

For r 7^ 0, the last term vanishes and, indeed, E is invariant to the substitution 
73(2) _^ Extending this field to the entire space, this last term contributes with a 

delta-singularity such that we must write 

E^^\r) = --^e,V,Ad,du- + ^kV5{r) . (32) 



Inserting equation (l2Ql) in equation (l32l) . and retaining only the delta-singularities, we 
can write 

E^g (r) = e,P,,9,5(r) + -^A V5(r) . (33) 

Substituting equation (l30l) in equation ( 1331) . one obtains equation ( 1291) . 

The advantage of employing the STF moments 'P^'^^ instead of the primitive 
moments P^"-* is manifest for higher order terms from the multipolar expansion. Even 
from n = 3, the contraction of the angular average of a product of more then six factors 
V with a primitive moment represented by a tensor which is only symmetric becomes 
cumbersome. 

Though, maybe, only of theoretical interest, let us consider the general case of 
arbitrary n. The STF projection 7"(7'*-"^) of the symmetric tensor P*^") is defined, up 
to a numerical factor, by the equation 

'T~ii...i„{^^ "*) = 'Pii...in ~ Pji...'t„ ~ ^{iii2^ii...in} i (34) 

where A^""^) is a symmetric tensor and is defined by the conditions of the trace-free 
character of T'^"). For low values of n (the ones of practical interest), the components 
Aii...in can be calculated directly from the equation system representing the vanishing of 
all the traces of the tensor 7-*'-"^ For higher orders n, there is a general formula known 
in literature [8], [12] which, with the notation from the present paper, is written as 

j(n:m) 



n/2] t_i^m(^2n - 1 - 2m)!! 

(2n- 1)!! ' 

^i2^i in denotes the components of the {n — 2 m)-th order tensor obtained from P*^"^ 



{n:m) 



m=0 



by contracting m pairs of symbols i. This equation is known as the detracer theorem 
[I2] . As a consequence of this theorem, the components of the tensor A^""^^ are written 
as 



[n,/2-l] 



-l)-[2n-l-2(m + l)]!!^ ^ 



A. . rp(n)l _ V -J i^x^ - ^yno^.,y. p(n:m+l) .ggv 

l\^...ln-2[^ J Z_^^ [m + l){2n - "1*1*2 ■■■ "«2m-l«2m'^i2m+l...jn-2}- ^"^^ ' 
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Inserting equation in the expression of E^^\r) given by equation fl24p . we obtain 

AireoW. r 47reo'T-' r 

and, since the contraction with the fully symmetric tensor A*^""^) produces = 
n{n — l)/2 identical terms, 

^ ^ 47reon! " r 87reo(n - 2)! ' r 

7?W||V"+^-+ / ' A""'||V""^(5(r) ■ (37) 



Searching the extension of to the entire space, for establishing the 5-type 

singularities, we have to calculate the limit 

0) = izl)!l^iiin / d^x fpW||V"+i-) 0(r) 
V (0)' ^/ 47r£nn! ^^o y^-^ V r/ ^ ^ 



-0 

+ 



2(n-2)!£o 



lim / d^a;A("-2)||V"~^<5(r)(/)(r) 



or 



' ^/ 47reo n! e^o Jv, \ ■ r 
1 



2{n-2)\eQ "V ^/o 

As done in Ref . P , the limit of the remaining integral from equation ( 138|) can be easily 
expressed for arbitrary n. Applying the Gauss theorem, we can write 



lim^ d\ (^pW||V''+^i) 0(r) 



lim e, 



.(39) 



Introducing the Taylor series for the function </)(r) in the surface integral from the last 
equation and since r = e on the sphere S^, we can write 

liniCi / d^ Vi (Vi,...i,A, . . . 0(r) 

= Att Ei lim ^ , ( C'g^„z/i„^, . . . (^i^^, . . . di^^^ (f))^ . (40) 



Let us evaluate the tensorial contraction which is present of the general term in the 
series from the previous equation: 

From equation ([7]) we can easily see that all the terms containing at least a symbol 5i^i^ 
with 1 < q, s < n give null results by contraction with the traceless tensor Only 
the term corresponding to /c = in equation ([7]) can give results different from zero for 
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this contraction. Moreover, a result different from zero can be obtained if and only if 
a + 1 > n. Therefore, equation fH2l) can be written as 



lime, / dSu, fv,,...^,A^---di^-) 0(r) = 47r e,(-ir(2n - 1)!! 



^ a=n-l 



X 

For a > n — 1, the corresponding terms from the series in equation f HTl) contain positive 
powers of e. Consequently, the corresponding limits for e — )■ vanish. From this series 
only the term for which 

a = n-l (42) 

can be different from zero and the result of the limit in equation f l4ip is given by 

1\ 47r(-l)"(2n - 1)!! 



lini Bi f dSui ( Vi,,„i„di^ ...d, 



rJ (n — 1)\ 



^ 47r(-l)-(2n-l)!! / ^(n)n2n\ II 

{n-l)\ \ " /"V^ 
Let us evaluate now the volume integral in equation (139|) : 

/■ / 7?("-) \ 
limy d?x (^"11^— J ^<P{r) 

dSi^W V'^Ml—- V(f)(r) - d^x V'^Ml—- V^(h(r) 



(43) 



lim i dSu\\[V^'\\- )V0(r)-/ d^x(V"~^||^ 1 V'0(r) .(44) 

£->o [Jse \ r J Jve \ r ) J 

The part corresponding to the surface integral can be written as 

lim / d5i.ll fv"-i||^^ V0(r) = e,P,,..,„ lim / d5i/,„ f^,, . . . 9,0(r) . 

(45) 

Introducing the Taylor series for 0(r) and standing out the average over i/, we obtain 

j- ( 73(")\ 

limj^^ dSv\\ ( V"-^||^j V0(r) 

= 47r lim —Vi^...i^ ( Cf~l^~^\i^ . . . u-,^^^) x . . . d,^^^di<p)^ . (46) 

Analogously to the reasoning from the previous case, we can see that the limit is zero 
since for a = n — 2 



il...in \ I'n ^n + l ''2n — 2/ 



Therefore, 



r / 'pH\ I- ( 'p(")\ 

lim / d^x V\\- — V0(r) = - lim / d^x — V20(r) , (47) 
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and by repeatedly applying the procedure, all terms cancel. Finally, only the surface 
integral from equation fl39|) gives a limit different from zero and 

hm / d'x fpHllV^+i-) 0(r) ^ 47r(-ir(2n-l)!! 2n\|| f^n-iA ^^g) 

e-^oJv, V rj ^ ' (n-1)! \ n / n V ^/o ^ ^ 

Let us consider the contraction 

To this contraction contribute only the terms from the average of the //-product not 
containing factors (5jj.jp with 1 < k, I < n. According to equation f|T5|l . the terms giving 
non-zero contributions are of the form 
1 



^iljl • • • ^injn 



{2n + l 

and all such terms are obtained considering the n\ permutations of the indices ji . . .jn 
in this product. Therefore, the final expression in equation fH5]) is given by 

hm / d^x fpWllV^+i-) 0(r) = l!^zl)^pW|| fv^-V) . (49) 
£^oJv, V rj ^ ' 2n + l ^ ^ 

Equation (l38l) . with equation ( H9l) inserted in, yields 



^!S,0) = -7 TTT^ — P(")||fv"-V) - , \. A("-^)||(V"-V) .(50) 

(°)'^/ (n-l)!(2n + l)5o ^ 2{n-2)\eo "v ^^o ^ ^ 

Therefore, the delta-form distribution associated to the 2"-polar electric field is given 
by 



= i-il! ^p{«)||v"-i5(r-) + _Lll-_A("-2||V"-i5(r-) . (51) 

(n - 1)! (2n + l)eo 2(n-2)!£o 



One can easily see that the result fl32|) for = 2 is, indeed, a particular case of the 
formula ( 1511) . In the case n = 3, equation ( 15T1) becomes: 

The result can be expressed in terms of primitive tensors P*-"^ introducing equation (E 
in equation f lSTj) : 



-1)" . (-1)" 

'2n- 



E^S. = 7 TTi^ pi'-)\\v^-H(r) + / A("-^||V"-^(5(r) 

(n- l)!(2n + l)eo 2(n-2)!eo 



-ej(5{jj,Aj2...j„_i}9,, . . . 9j„_i5(r) . 



(n- l)!(2n + l)£o 

In the last term from this equation, there are n — 1 identical terms of the form 
and C^_i = {n — l){n — 2)/2 terms of the form 
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Finally, 

- (n-l)!(2n + l)eo " ^ + 2(n - 2)!(2n + l)eo " ^ ^ 

— 7 TT7^ — A(""2)||v"-3^(5(r) . (52) 

2{n - 3)\{2n + l)eo n V ; v ; 

This last equation becomes equation (!29|) in case n = 2. The terms containing A5(r) 
are present beginning from n = 3. 



4. Singularities of the magnetostatic field 

For the vector potential in the exterior of the domain V, we have 
A(r) = V V X fv"-i||M(") 

= E ^P^n • • • d..-."^^"^, (53) 

where M*^*^) is the magnetic n-th order moment defined by the Cartesian components 



ft C 

Mn...i„(t) = ^^y^d^xx,, ...a;i„_,(r X J(r,t)),^, (54) 
or, with tensorial notation: 

M(")(t) = j d^x X J(r,t) . 

The corresponding expansion of the magnetic field B{r) = V x A{r) is given by 

B(.)^i^vxE^vxfv-||M!:l) 

47r n! ^ W 



E 



47r ^ n 

n>l 



V- V"|| - A V 



(55) 



For r 7^ 0, the last term containing A(l/r) is not contributing to the multipole expansion 
since A(l/r) = 0, but, searching the extension of this expansion to the entire space, 
including the point O, we have to consider it. This term, extended as in the electrostatic 
case, as An M<^")||V" ^5(r), is considered as a first extension of B. It remains to process 
the limit that implies the expression V (^M*-"''|| V"(l/r)^. 
In the dipolar case, we write 

S(i) = if^ei9,a^ + 47rm5(r) 

and, applying equation f lT6l) . we obtain for the delta-singularity, the well-known 
expression [2]: 

"^(0) ^ ~^"^'^(^) + /^o"^^(^) = ■^^m5(r) . (56) 
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For the higher multipolar orders, let us begin with the 4-polar term: 

8n \ r J Svr r 2 

Applying equation fl20|) . 

since SjkMjk = Mjj = 0. Further, 

Sgj = ^e. M,, d,6{r) - ^eM,^d,5{r) = ^ m'^ \\V5{r) - ^ (V5(r)) ||M(^) , (57) 



^ (2) 

where M is the symmetric part of the tensor M*^ corresponding to the identity 

= ^ + Mji) + ^ (M,, - Mji) =Mij +^ei,kMk ■ (58) 

The antisymmetric part is expressed in terms of the components of a first rank tensor: 

2 



Ni = EijkMjk = - I d^a; [r X (r X J)]. . 
6 Jv 



For n = 2, the symmetric projection of the magnetic moment is an STF tensor, i.e, 
denoting by Al^") the STF magnetic moments, 

the equation fl58l) is written as 

1 

The introduction of these result in equation (157|) gives 

^(0) = -^-^M^''\\VS{r) - ^AT X V5(r) , (60) 
where N = NjCj. 



M,j- = A^i^- + -Eijk^k ■ (59) 



Beginning from n = 3, the symmetric part of the magnetic moment tensor is not 
the same with the STF one. We write the identity 

where Mj^j^is is symmetric in the first two indices. The first parenthesis represents the 
fully symmetric part of the tensor M'-'^^ and the second one can be expressed in terms 
of the second order tensor N^^) defined by the components 

with the relationships: 



We can write 



^Mj-^jjis ^i2i3h) ~^ ('^ii«2«3 ^hhh) ~ ^hisQ^hq ~^ ^i2i3q^iiq ~ 5^ ^i\i-iq ^{hi2q) ' 



A=l 
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where by the notation ^[i^i^q) we understand the component without the ir 
Employing the above definitions and notation, equation ( l6Ti) can be written as 

1 
3 



where by the notation ^[i^i^q) we understand the component without the index ix. 

LUC 

^ 1 ^ (A) 



A=l 

^ (3) 

where M is the symmetric part of the tensor M*^ ^ 

Let us write the third order term from the magnetic field expansion (l55l) 

^ ' 247r V r J 247r " V r J 

Retaining the second expression which represent an extension as distribution with the 
point-like support O: 

B'"Hr) = (v^'ll^) + ^ {V%r)) ||M<') , (63) 

^(3) 

Introducing the symmetric tensor M from equation f l62|) in the first expression of the 
right-hand side of the above equation, 

/ ^ (3) \ 

Since in the last expression, the two terms of the sum contain either the di^di^Si-^i^q or 
di^di^Si^i^q which vanish, we can write 

/ (3) ~ 

V|V»||M^!l).vMl^^^| . (64) 
This result expresses the invariance of the multipole expansion of the magnetic field to 

^ (3) ^ (3) 

the substitution M M . The introduction of M in the 5-type singularity from 
equation (l63l) gives 

^ (3) 1 2 
A=l 

X «(3) 2 ^(3) „ 2 

V'5(r)j II M +-eiHi^qeiqiAiAi^5{r) =M || V'5(r) + -eini^qSiqi^Ai^Sir) 

•H- 

since M is symmetric. For expressing in a compact form such tensorial contraction, let 
us introduce the notation 

^1 X |B^ ^ ^i/^ii...in-iq^iqs^ s ii...in—\ ■ (6^) 

Then, 



(V2(5(r)) ||M(3) =M ||V25(r) + -N(2)| X |V'(5(r-) . (66) 



^ ^ 247r 
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The octupolar term B^^'' including partially 5-form singularities becomes 

/ o (3) \ 

'=^11^ + ^ M^'^ \\V'6ir) + ^N(2)| X iV'Sir) . (67) 
y r j 6 9 

The STF projection up to a numerical factor, is given by 

« ^ ^ (3) 

Mi^i2i:, =Mj^i2i3 -5{jjj2Aj3}(M ) , (68) 

where the symmetric tensor A*^ ^ corresponds to M''"^ by the a formula of the type 
(151|1 . It easy to see that 

~ 1 1 1 /" 

Ai = - Migg= —Mgg, = — / d^x (r X J)i . (69) 
The introduction of equation fl68l) in equation fl671) gives 

i^4r7r \ y o o 

+ ^N(2)| X |V25(r) . (70) 
9 

From this last equation, it is seen that the multipole expansion of B{r) (r ^ 0) is 
invariant to the substitution M^^^ Al^^^ P fTO], 

It remains to calculate the extension of the first term from equation (170|) to the 
entire space. For this, we have the result ( H9l) from the case of electrostatic field which, 
for arbitrary n, gives in the present case: 

hm / d^xv(v-\\^] 0(r) = ^"(-^)"" a4(")|| (V"-V) • (71) 

The final result for the singular part of B^^^ having as support the point O is: 



B 



(3) 



_ ^>t(=')||v25(r) - l^A\\V^6{r) + ^A A5(r) + ^N^^)] x \V^6ir) . (72) 
2 X 6 6 9 

Let us consider the extension for an arbitrary n of 

BW(r) = f Vll^l + /f£izi)!:i (V""^5(r)) ||MW . (73) 



The generalized equation (162|) is given by 



-j^ 71—1 



M,....„ =M,...,„ +^ 5: e.,.„,N[f) , (74) 



„ . . - (n...in-l'?) 

A=l 



where it is introduced the (n — l)-th order tensor N^" ^\ partial symmetric in the first 
n — 2 indices and with null contraction of the last index with any of the indices 

iq, q < n - I: 



^ f 3 

^i„_ipsMj^,,,j^_2Ps = ^ ^ ^ j^d xxi^. . .Xj^_2 [r X (r X «/)]j^_^ • (75) 
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Inserting equation fl74p in the expression of B^^^r), we consider the different terms 
from equation fl73l) : 



\ 



A=l 



(n) ' 

M 



since all the terms from the last sum contain a contraction of the type £iii„qd^di^ with 
1 < / < n - 1. 

n-l 



n 



A=l 



/ 1 \ ri — 1 

(V" '5(r)j II M + ei ni^...i„^,q£iqiAA2 ■ • • ^n-i5(r) 



n 

n — 1 , 



M \\V'^6{r) + N("-i'| X |V"-^5(r) 



using the notation ( 165|) . The final result for the substitution of equation (174|) in equation 
T3l) can be written as 



/"o(-l) 



n-l 



-V 



M 



+ 



Ann). 

/io(-l)"^^ (n-l) 
n! n 



n\ 



M ||V"~'(5(r 



(76) 



Introducing the STF tensor Al*-"-' with the components given by an equation of the 
type dMl) 



M 



and considering the different terms from equation (1761) . we obtain 



V 



V V 



„ A<(")\ n(n-l) 



>13...2„ 



V V 



27rn(n- 1) A^" ^\\V''-^6{r) . 



The second tensorial contraction from equation (|76l) can be written as 



M ||V"-i5(r) = A4(")||V"-i5(r) + e,„9.,...9.„_,5(r)%,,A,,.,„} . 

In the last expression, there are n — 1 terms containing the factor 5jqi„ , q = 1, ... n — 1 
and C^_i = (n — l)(n — 2)/2 terms containing the factor 6i^i^ with q and s between 1 
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and n — 1, such that the last equation can be written as 

■in) 



+ 



{n -l){n- 2) ~(n-2 



V"-^A(5(r) . 



(77) 



Collecting all the above results, equation ( l76l) can be written as 



A V"+^ 1 1 !_!_ + 1 1 V"-^(5(r) 
47r r 



- (^-l)(^-2)^(n-2)^j^„,,^^^^ ^ (n-l)(n-2)~(n-2)^^^(„_3)^^^^^ 

n — 1 



+ 



n 



(78) 



Employing equation fl7T]) . we finally write the singular (5-form part of B^'^'^ 



S(")(r) 



+ 



2n + 1 

(n- l)(n-2) ~(n-2),,^(„ 



7T. ..,'„-) II _„„i ^ ("■ —!)("' — 2) ~(n-2) -. 



V("-=^) A5(r) + - — 1n("-i)| X |V"-^5(r) 



(79) 



2 ■ ' n 

One can easily verify that equations ( l60l) and (1721) are particular cases of equation ( 1791) 



5. Conclusion 

The results of the present paper concerning the (5-form singularities of the 
electromagnetic field in the static cases are not so appealing as the ones done in Ref. 
[1]. However, in [1] the corresponding results are obtained, in our opinion, without 
employing the full content of the multipole expansions. The unpleasant presence of the 
parameters A and N in the expressions of the 5-type singularities is a consequence of the 
hypothesis that the basic multipolar moments are the primitive ones (p(") in the electric 
case and M^") in the magnetic one). The employment of irreducible representations 
by the STF tensors has the advantage of simplicity in expressing some quantities in 
several circumstances. It seems, for example, that some trouble appears when employing 
the multipole expansions in spherical coordinates. Let us consider the example of the 
electrostatic potential multipole expansion: 

oo I 



Mr) 



1 



Q 



47re 



l=Om=~l 



and the particular case 1 = 2: 



E '52m ^2 



, •3 ^"i - 2m ■ (80) 

The spherical moments Q2m are linear combinations of the components Vij of the STF 
moment 'P^'^\ Writing this term in Cartesian coordinates, 

$(2)(r) 



^Lp(2)||v2l 

87r£:o r 
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and introducing the STF tensor V^'^\ 



Stt^o 
i.e 



Viididi- + AA- 



For the first term, equation ( 1161) gives the 5-singularity 

since Vu = 0. Therefore, the electric dipolar potential has a delta-type singularity 

$(2)) = — -A6(r). (81) 
ho) 2^0 

But, for the term expressed in spherical coordinates (equation (IHOj) ). since there are no 
derivatives of 1/r, one has no (5-type singularities. 

In our opinion, equation flHTj) represents the correct result. From the given example 
we also see that some care is necessary when the invariance properties of the multipole 
expansions are used. 

If correct, the results of the present paper can be useful in classical and quantum 
physics, in the second case starting with the problem of the hyperfine atomic structure. 
The dynamic case will be treated similarly elsewhere. 
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